On magnetic field generation in Kolmogorov turbulence 
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We analyze the initial, kinematic stage of magnetic field evolution in an isotropic and homoge- 
neous turbulent conducting fluid with a rough velocity field, v(l) ~ l a , a < 1. We propose that 
in the limit of small magnetic Prandtl number, i.e. when ohmic resistivity is much larger than 
viscosity, the smaller the roughness exponent a, the larger the magnetic Reynolds number that is 
needed to excite magnetic fluctuations. This implies that numerical or experimental investigations 
of magnetohydrodynamic turbulence with small Prandtl numbers need to achieve extremely high 
resolution in order to describe magnetic phenomena adequately. 



1. Introduction. In a turbulent highly conducting 
fluid, magnetic fields may be amplified since the field 
lines are generally stretched by randomly moving fluid 
elements in which these lines are frozen [1] . Such a mech- 
anism of turbulent dynamo is expected to work in a vari- 
ety of astrophysical systems (galaxy clusters, interstellar 
medium, stars, planets), is confirmed numerically, and is 
consistent with simple analytical models. 

Valuable insight can be gained from considering the 
so-called kinematic stage of the dynamo, when magnetic 
field is amplified from an initially weak 'seed' field. As 
the simplest example, assume that the resistive scale is 
smaller than or of the same order as the viscous scale 
of the fluid. In the turbulent Kolmogorov velocity field, 
the smallest eddies have the highest shearing rate given 
by v(l)/l ~ l~ 2 l z , where I is the size of the eddy, and 
v(l) ~ I 1 / 3 is its turn-over velocity. Therefore, at this 
stage the magnetic field grows predominantly on small 
scales. Fig. 1 shows schematically the initial stage of the 
evolution of the magnetic spectrum, see, e.g., [2]. 

The magnetic energy collapses toward small, resis- 
tive scale during this initial evolution, until the field is 
strong enough to affect the dynamics of fluid through 
the Lorentz force. Such behavior is the evidence of 
small-scale turbulent dynamo; it is firmly established 
in numerical experiments, and can be derived analyti- 
cally [2,3,5-9]. Since in this example, the resistive scale 
is smaller than the viscous one, the dynamo is essentially 
governed by a smooth, viscous-scale velocity field. 

To study the growth of magnetic energy on larger 
scales we have to understand how the magnetic field is 
generated in the inertial interval of the turbulence (the 
interval of scales much smaller than the external scale 
where the turbulence is excited, and much larger than the 
viscous scale where the turbulent energy is dissipated). 
This question is non-trivial since in this interval the ve- 
locity field is not smooth, i.e., v(l) ~ l a , with a < 1. 
This situation is especially relevant for the case of small 
magnetic Prandtl number (Pm = v/r), where v is fluid 
viscosity and r\ is resistivity), where the magnetic energy 
is concentrated mostly in the inertial interval of the ve- 
locity field. 

The question was first addressed by Batchclor [3] who 



used the analogy between magnetic field lines and fluid 
vorticity lines to conclude that when v < rj, magnetic en- 
ergy is not amplified. This analogy was criticized in [4] 
because the magnetic field can have arbitrary initial con- 
ditions, while the vorticity field is related to the veloc- 
ity field. It was further argued in [10] that magnetic 
field line stretching generally dominates resistive recon- 
nection, thus making dynamo possible for v < r\. Direct 
numerical simulations of MHD turbulence with Pm = 1 
confirm that week magnetic fluctuations are generally 
amplified by the Navier-Stokes velocity field except for 
a special case when the initial magnetic field is close to 
the vorticity field [11,12]. 




FIG. 1. Sketch of the initial stage of magnetic field am- 
plification by turbulent velocity field. Ey. is the spectrum of 
the velocity field in k space, Hk is the spectrum of the mag- 
netic field. The viscous cut-off k v is larger than the resistive 
cut-off k v . Magnetic field energy is shown for three consec- 
utive times, it grows and gets accumulated at the resistive 
scale. 

Recently, there appeared the number of numerical sim- 
ulations of MHD turbulence with small magnetic Prandtl 
numbers, where magnetic fluctuations were not ampli- 
fied [13-16], which revived the claims that dynamo does 
not exist in the Kolmogorov turbulence with Pm <C 1, 
see e.g., [15,16]. However, astrophysical observations 
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show that magnetic fields are generated by turbulent 
motion rather effectively in planets and stars where 
magnetic Prandtl numbers are small (e.g., in the geo- 
dynamo, Pm <~ 10~ 5 — 1CP 6 , in the solar photosphere, 
Pm <~ 10~ 7 ). This apparent contradiction motivated our 
interest to the problem. 

In this paper we argue that dynamo always exists in 
a rough velocity field, such as the Kolmogorov field. 
We however find that the magnetic Reynolds number, 
Rm = v(Lo)Lo/r], and the numerical resolution required 
to generate magnetic field strongly depend on the rough- 
ness exponent of the velocity field, a; the rougher the 
velocity field, the larger the required resolution. This 
result explains why dynamo action is hard to achieve in 
experiments with small magnetic Prandtl numbers, while 
it is easily achieved when the magnetic Prandtl number 
is large, Pm > 1. In the latter case the velocity field is 
effectively smooth at the resistive scales, where the mag- 
netic energy is concentrated. 

Our analysis is based on the so-called Kazantsev model 
of kinematic dynamo, which is exactly solvable and which 
allows us to change the velocity roughness exponent, a. 

2. The Kazantsev model for a rough velocity field. 
Consider the induction equation for the magnetic field: 

dB l /dt + v j dB i /dx j - B j dv i /dx j = V AB\ (1) 

where B l (x, t) is the magnetic field, v l (x, t) is the fluid 
velocity, rj is small ohmic resistivity, and we sum over the 
repeated vector indices. At the initial stage of evolution, 
the weak magnetic field is passively advected by the fluid. 
Therefore, one can use simplifying assumptions about the 
fluid velocity to make the problem manageable. Kazant- 
sev [5] and Kraichnan [17] introduced the model based 
on the Gaussian, short-time correlated velocity field, with 
zero mean and the covariance 

(v < (x,i)w 3 '(x',t / )) = ^'(x-x')<K*-*')- ( 2 ) 

This model is a valuable tool for the analytical inves- 
tigation of kinematic dynamo; direct numerical simula- 
tions reveal that a purely Gaussian velocity field ampli- 
fies small magnetic fluctuations in a similar manner as 
the true Navier-Stokes field [11]. 

Assuming isotropy and homogeneity, the velocity cor- 
relation function has the form 

513 -— j+^M— . ( 3 ) 

where r = x x'. If we further assume that the velocity 
field is incompressible, we have kn — Kl + and 
velocity statistics can be characterized by only one scalar 
function, Kl{t). 

The model defined by (1), (2), and (3) allows one to 
write a closed equation for the correlation function of the 
magnetic field, 

(B i (x,^(x',i))= J ff«(x-x',i), (4) 



where, analogously to (3), the H l i function can be rep- 
resented as 

S^-—)+H L (r,t) — , (5) 

furthermore, the condition V • B = gives Hjy — 
Hl + (rH' L )/2. We will characterize the magnetic field 
correlator by the function Hl {r, t) . The equation for this 
function can be found by differentiating (4) with respect 
to time, and by using Eqs. (l)-(3). A rather tedious but 
essentially straightforward calculation gives: 

d t H L = kH'I + (^-k + n^j H' L + (k" + ^ H L , (6) 

where primes denote the derivatives with respect to r and 
we have introduced the 'renormalized' velocity correla- 
tion function n{r) = 2ti+kl{0)~ kl(v). Equation (6) was 
originally derived by Kazantsev [5], and can be rewrit- 
ten in a different form that formally coincides with the 
Shrodingcr equation in the imaginary time. Effecting the 
change of variable, Hl = ip(r, i)r _2 /t(r) -1 / 2 , one obtains: 

dtip = K(r)ip" - V(r)ij>, (7) 

which describes the wave function of a quantum parti- 
cle with the variable mass, m(r) = 1/[2k(t)], in a one- 
dimensional potential (r > 0): 

V(r) = £«(r) - \n"{r) - ^'(r) - (8) 

This equation can be investigated for different choices 
of K,(r), however, we restrict ourselves to the inertial in- 
terval of the turbulence, where the velocity correlator 
has power-law asymptotics, n{r) cx . The exponent 
[5 can be found from the scaling of turbulent diffusivity, 
D r~j v(r)r <~ r 1+a . Indeed, in the derivation of Eq. (6) 
we used the integral D = f£°([v(x,t) — v(x' , t)][v(x, t') — 
v(x' , t')])d(t — t') = n{r), which is the turbulent diffu- 
sivity [10]. Comparing the two expressions we obtain 
0=l + a. 

The Shrodinger equation (7) has the effective poten- 
tial U eff (r) = V(r)/K(r) = A((3)/r 2 , where A((3) = 
2 - 3/3/2 - 3/3 2 /4. When A([3) < -1/4, the quantum 
particle falls toward the origin [18] and its wave function 
is therefore concentrated at the smallest, resistive scale. 
This behavior is the manifestation of the dynamo mech- 
anism that we discussed in the introduction. We obtain 
that A(l + a) < —1/4 for any roughness exponent of the 
velocity field, < a < 1, therefore, dynamo is always 
possible; the same conclusion was made in [10]. 

3. Magnetic field correlator and dynamo growth rates. 
To find the wave function (the magnetic field correlator) 
we need to specify the boundary conditions. The small- 
scale regularization is naturally given by ohmic resistiv- 
ity. For scales much smaller than the correlation scale 
of the velocity field, the k function can be expanded as 
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/■c(r) ~ kl(0)(2t7 + r 1+Q ), which corresponds to the limit 
of infinitely small Prandtl number. In this formula and 
below we use the dimensionless variables: r\ is measured 
in the units of the large-scale turbulent diffusivity, kl(0), 
while r is in the units of the integral scale, L . 

The boundary condition at the origin follows from 
fmitcness of magnetic energy, lim Hl(t, t) = Holt) < oo. 

r—*0 

The boundary condition at large scales follows from the 
absence of the mean magnetic field, lim Hi,{r,t) = 0. 

r — >oo 

We will see in a moment that at the kinematic stage mag- 
netic energy is concentrated at the resistive scale, and the 
corresponding wave function decays exponentially fast in 
the inertial interval, r 3> r v , so it is indeed independent 
of the large-scale properties of the velocity field, as we 
expected. 

Problem (6) can be cast into the Sturm-Liouville form 
by changing H^r.t) = h(r,t)/r 2 , which teaches us that 
its solution can be expanded in the cigenfunctions of the 
corresponding Sturm-Liouville operator. The maximum 
growth rate is therefore given by the largest eigenvalue of 
the operator on the right-hand side of Eqs. (6) or (7). The 
corresponding wave function is the ground state of the 
potential U e ff(r). Thus, following Kazantsev [5], we look 
for the solution of (7) in the form tp(r, t) — <p(r) exp (At). 

In the inertial interval, 2r) <C r 13 , the equation for the 
ip function reads 



-<p + 



-3/3 2 -6/3 + 8 X_ 
4r 2 r@ 



<P = 0, 



(9) 



where A is dimensionless and is measured in the units of 
kl(0)/Lq. For small A, the potential is dominated by the 
first term in the square brackets, and one can easily check 
that the corresponding wave function oscillates in the in- 
ertial interval. However, the ground-state wave function 
cannot oscillate, therefore the ground-state growth rate, 
A, will be such that the second term in the square brack- 
ets of (9) dominate the first one in the inertial interval. 
At the resistive scale, r n = (2?y) 1 /' 3 , both terms should 
be of the same order, therefore, A ~ r/il 3 - 2 )/! 3 . 

The wave function corresponding to the growing so- 
lution, A > 0, decays exponentially fast for r 3> T] 1 ^, 

ip oc exp ^— 2v/ ^2-^ g)/2 ) • ^his mnc tion is concentrated 
at the resistive scale and its growth rate is of the order of 
the eddy turn-over time at this scale; in agreement with 
our qualitative discussion in the introduction. When one 
changes the magnetic Prandtl number (by changing vis- 
cosity, for instance) the effective roughness of the velocity 
field at the resistive scale changes. We therefore suggest 
that the effect of different Prandtl numbers can be stud- 
ied in terms of Eq. (9) with different roughness expo- 
nents. By this analogy, the limit of the smooth velocity 
field, (3 — 2, corresponds to large Pm, while the other 
limit of the Kolmogorov-scaled velocity field, [3 = 4/3, 
corresponds to small Pm. 

We observe that the rougher the velocity field, the 
broader the wave function compared to the resistive scale, 



and, therefore, the larger the magnetic Reynolds number 
necessary to generate magnetic fluctuations. To estimate 
the critical resolution for the dynamo onset, we solved 
Eq. (6) numerically with the large-scale boundary con- 
dition, Hl{1) = 0. For given (3 and i], we increased the 
'system size', £, until Hl started to grow and we thus 
found the critical value I = £ c . An analogous calculation 
could be done by fixing t and changing 77. 

To characterize the inertial range, we introduce the 
resolution parameter, R = l c jr n = £ c /(2r]) 1 ^; this pa- 
rameter is universal, i.e., it is independent of large and 
small-scale properties of the velocity field. We obtain 
that when the velocity roughness increases from (3 = 2.0 
to (3 = 1.3, the corresponding resolution parameter in- 
creases from R w 3.8 to R w 29, see Fig. (2). This 
means that if we simulated the random equation (1) di- 
rectly, the required numerical resolution should increase 
by about an order of magnitude in each spatial direc- 
tion as we go from the smooth to the Kolmogorov-scaled 
velocity field. 




1.7 1.6 1.5 

ROUGHNESS EXPONENT, BETA 



FIG. 2. The numerically found resolution parameter, 
R — £c/r v , as a function of the velocity roughness expo- 
nent, (3. 

The subresistive scales, r -C r/ 1 ^ , can be easily de- 
scribed using Eq. (6). By direct substitution, one can 
check that the magnetic field correlator can be expanded 
as H L (r,t) = H [l - r^/(2r]) + ...]exp(Ai), which im- 
plies that the spectrum of the magnetic field decays 
as Hk ~ k 2 \Bk\ 2 cx fc~ 1 ~' 3 in the subresistive region, 
Ty- 1 // 3 <§; fc <C v~ x IP . With such a spectrum, the rate 
of magnetic energy dissipation at the viscous scale is 
bigger than the rate at the resistive scale by the factor 
kl\B kv \ 2 /kl\B kv \ 2 ~ Pm- {2 - fi)/l3 . Thus, the effective re- 
sistivity may be 'anomalously large' in turbulent flows 
with small Prandtl numbers. 

Strictly speaking, the last result requires short-time 
correlated small-scale eddies and is not applicable in the 
Kolmogorov turbulence where velocity is correlated at 
the eddy turn-over time, T corr (l) <~ l/v(l). At the sub- 
resistive scales, I < l v , this correlation time is larger 
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than the resistive relaxation time, t v ~ I 2 /rj, while in the 
Kazantsev model the velocity correlation time is smaller 
than the resistive relaxation time. However, model (2) 
is physically self-consistent. The large magnetic energy 
dissipation is balanced by the large energy transfer from 
the small-scale velocity eddies to the small scale magnetic 
field; this is possible since the fluid energy is formally 
infinite in (2). In a practical situation where the veloc- 
ity correlation time is not infinitely small, the asymp- 
totic tail Hk ~ fc -1- ^ will hold up to the scales where 
this correlation time is comparable to the resistive time, 
k ~ l/^/?7T corr . One can expect that for smaller scales, 
the magnetic energy spectrum will have a steeper decay, 
see, e.g, [4,19,20]. 

4. Conclusions. We proposed that magnetic fluctua- 
tions are always generated in isotropic and homogeneous 
three-dimensional turbulence if the magnetic Reynolds 
number is large enough. The required critical magnetic 
Reynolds number sharply increases with velocity rough- 
ness, which explains the 'no-dynamo' outcomes reported 
in numerical simulations with small Prandtl numbers, 
see e.g., [13-16]. We do not see any physical reason 
for the absence of dynamo in these numerical simula- 
tions rather than lack of resolution. Our results also sug- 
gest that obtaining small-scale dynamo will be a serious 
challenge for laboratory experiments, where the magnetic 
Prandtl number is small, Pm ~ 10~ 5 , while the magnetic 
Reynolds number is rather moderate, Rm ~ 100 [21-24]. 

As we mentioned in the introduction, the Batchclor 
analogy of magnetic field lines and vorticity lines would 
hold for v = n if the initial magnetic field were propor- 
tional to the vorticity field. In this special case, mag- 
netic energy would not be amplified. In the model we 
investigate there are infinitely many such special initial 
conditions; any magnetic correlator Hr,(r) lacking the 
growing cigenfunctions in its expansion will not be am- 
plified. This, of course, does not mean that dynamo does 
not exist. 

We also note that the Schrddinger-type equation of the 
form (7) was considered for the magnetic field correlator 
in the case of a general velocity field in [10]. In this sense 
the applicability of the equation may not be restricted 
to the Gaussian, short-time correlated velocity field only. 
For example, we also investigated the modified Kazantsev 
model, with a more realistic, finite-time correlated veloc- 
ity field. The velocity correlation time at a given wave 
number, k, was chosen to be of the order of the eddy 
turnover time at the corresponding scale. This model 
could not be solved analytically, however, its numerical 
integration gave qualitatively the same results as (1,2). 
We will report these results elsewhere. 
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